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This paper studies solute segregation and its relationship to grain boundary energy in binary alloys. Grain 
boundaries are simulated using a binary alloy phase field crystal model that incorporates atomic elasto-plastic 
effects on diffusional time scales. Grain boundary energy versus misorientation data is semi-empirically fit 
to Read-Shockely theory in order to characterize its dependence on the average alloy concentration and un- 
dercooling below the solidus. The Gibbs Adsorption Theorem is then used to derive a semi-analytic function 
describing solute segregation to grain boundaries. Our results are in good agreement with direct simulations of 
the binary alloy PFC model. We also investigate how size mismatch between different species and their interac- 
tion strength affects segregation to the grain boundary. We interpret the implications of our simulations on the 
material properties related to interface segregation. 



I. INTRODUCTION 

Microstructure in metals is important for determining 
many of their properties (e.g., mechanical, thermody- 
namic, electrical). The various defects associated with 
microstructure formation (e.g., grain boundaries, dislo- 
cations, vacancies) contribute to an excess of free en- 
ergy of a system. As a material evolves towards equi- 
librium, its microstructure changes and along with it the 
material's properties. Grain boundaries are among the 
most important defects in metals. Their energy, com- 
position, and distribution directly affect the flow of dis- 
locations and influence the thermodynamics of second 
phase and precipitate formation. Being able to deter- 
mine grain boundary properties is therefore crucial to be 
able to properly model the myriad multi- scale properties 
linked to grain boundaries in solidification and solid state 
transformations. 

In an alloy, segregation of solute atoms can alter grain 
boundary energy CHS). The effect of segregation can 
also manifest itself in other ways. Two other properties 
strongly affected by solute segregation are solute drag 
[HO and grain boundary wetting El EHHl - In the for- 
mer case, the grain boundary energy is reduced by so- 
lute segregation, thus reducing the driving force to re- 
duce surface area (excess free energy) of a grain bound- 
ary. In the latter case, solute segregation can dramatically 
affect the thermodynamics of grain boundary formation; 
not only can segregation alter at what undercooling grain 
boundary wetting occurs, but it can allow for different 
grain boundary states (e.g., grain boundaries widths) (7). 
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There have been a number of experimental studies of 
grain boundary energy involving pure materials l9lfTT1l 
and alloys (HEU- Many studies have focused on char- 
acterizion of solute segregation and distribution (21 [3) as 
solute segregation typically has an important effect on 
grain boundary energy as demonstrated in (T] [TJl [131 . 
For pure materials and dilute alloys, many of these stud- 
ies have found that the grain boundary energy is well- 
fit by the well-known Read-Shockley Law when neigh- 
bouring grains are misoriented by small angles l9l fTT1l. 
It is also possible to adjust the parameters of the Read- 
Shockley equation to fit a larger range of misorientation 
angles I9UTT1. 

A number of theoretical approaches have been con- 
sidered to determine grain boundary energy in metals. 
The most prevalent, for both pure materials and alloyed 
metals, are the analytic analytic and semi-analytic dislo- 
cation models of Read and Shockley ['11, 14 ] and Van 
der Merwe Q, and models employing simple thermo- 
dynamic considerations of an interface [12) . Various 
computational approaches have also been employed to 
determine grain boundary energy in pure metals, includ- 
ing Monte Carlo simulations ITT31 and lattice statistics 
fT6) . Some computational approaches have also been 
used to model solutal effects in grain boundaries. These 
include monte carlo methods 0, molecular dynamics 
l3l[T3l and phase field simulations Q. 

Among the continuum theories available to model so- 
lidification and solid state transformations, the phase 
field crystal (PFC) method has recently emerged as a 
promising method to access the time scales and ef- 
ficiency of traditional phase field approaches while 
self-consistently incorporating many salient features of 
atomic- scale elasto-plasticity inherited form its connec- 
tions to fundamental microscopic theories [17] [18) . Per- 
haps the most important impact of the phase field crystal 
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approach will be its ability to bridge the divide that ex- 
ists between microscopic models and traditional phase 
field models lfT9l . For example, through the application 
of coarse graining theories, several researches have used 
PFC theories to compute the magnitude and and orienta- 
tional dependence of crystal-melt surface tension in pure 
materials (20|[22 and alloys l22l . PFC simulations of 
grain boundaries of pure materials have also convinc- 
ingly elucidated the Read Shockley behaviour for grain 
boundary energy versus misorientation (23l |24l and de- 
tails of pre-melting in grain boundaries I24ll25ll . 

In this work, we will use a recent phase field crystal 
model of binary alloys 1 17 ] to study the mesoscale be- 
haviour of how alloying and undercooling affect solute 
segregation and grain boundary energy in alloys. Nu- 
merical simulations in conjunction with an analysis us- 
ing a a single-mode decomposition of the crystal density 
field will be exploited to develop an expression for grain 
boundary segregation by exploiting a modification to the 
Read- Shockley equation that incorporates the effect of 
impurity concentration and temperature. We also con- 
sider the effects of lattice mismatch on segregation to the 
grain boundary and the Cottrell atmosphere around a dis- 
location core. 



II. PHASE FIELD CRYSTAL MODEL OF A BINARY 
ALLOY 

The phase field crystal (PFC) model of alloys used in 
this work is derived in detail in Ref. ifTTl . The resul- 
tant PFC free energy is expressed in terms of a tempo- 
rally coarse-grained normalized crystal density field and 
a relative density difference that is analogous to a solute 
concentration field. In particular, the normalized total 
density is given by n = (p — pi)/ pi and the normal- 
ized concentration by ip = (pi — P2)/ Pu where the to- 
tal density p is the sum of the density of each species, 
p — pi + p2, and pi is the density of a reference liquid 
state. The Helmholtz free energy functional expressed in 
these variables is given by 

F = ^{(^ + B 2 V)y +J B X n(2V 2 + V 4 )| 

+ rjB x nil)(\7 2 + V 4 )n d 3 f} (1) 

where Bq is the isothermal compressibility of the liquid 
at ip = 0, B2 determines how the isothermal compress- 
ibility of the liquid changes with ijj, B x is related to elas- 
tic constants in the solid, and t, v, u are determined by 
the low order terms of a local expansion of the classi- 
cal density functional theory description of the material, 
w is related to the various atomic bond energies, and K 
is related to w and the lattice spacing IfTTl . The differ- 



ence Bq — B x plays the role of normalized temperature 
variable. All lengths are scaled such that the lattice con- 
stant is a = 4-7r/ >/3 when the lattice mismatch parameter 
r] = 0. The lattice spacing changes with concentration 
according to the parameter 77 = (I /a) da /dip. 

Assuming conserved dissipative dynamics for both 
fields, the evolution equations are: 

In Eqs. [2] and |3j the atomic mobilities have been ab- 
sorbed in the time variable and a noise term reflecting the 
effect of thermal fluctuations on the evolution of the sys- 
tem has been neglected. The chemical potentials corre- 
sponding to each conserved field are given by p n =SF/ Sn 
and p^=SF/Sip. 

Equations [T][3] can be represented on mesoscales by 
a set of complex order parameter equations, the coeffi- 
cients of which are directly linked to those of the above 
PFC model, which is, in turn, linked to a simplified clas- 
sical density functional theory of freezing. The complex 
order parameter representation of Eqs.[T][3]has also been 
shown to reduce to the form of a traditional scalar phase 
field model with coupled strain effects fT9l . To the ac- 
curacy of a single-mode approximation, such an analysis 
thus yields a microscopic connection between continuum 
elastic effects and solute concentration and temperature. 



III. THERMODYNAMICS OF GRAIN BOUNDARIES 

To investigate the thermodynamics of segregation be- 
haviour at a grain boundary in a binary alloy, we need 
to consider grain boundary energy. Grain boundary en- 
ergy in pure materials can be determined by a number 
of different methods. For grain boundaries with a small 
misorientation angle, Read and Shockley derived the re- 
lation now named after them, 

lgb = E o 0(A-\n(0)), (4) 

by considering the grain boundary as an array of disloca- 
tions, where the dislocation cores do not overlap. For a 
2D crystal, the constants in Eq.|5]are Eq = Y 2 b/ (87ra), 
where Y2 is the 2D elastic modulus, a = \/3/2 is a cor- 
rection factor for hexagonal, as opposed to square, ge- 
ometry, b « a is the Burger's vector of the dislocation, 
andA = l + ln(a/r )-ln(27r) « 1.5 - In (2ir), which is 
related to the core energy by the core radius, ro I23ll24l . 
For high angles, some theoretical approaches consider 
the forces between atoms in a fixed geometry (21 [T6l . 
while others treat the high angle grain boundary as an 
amphorous phase sandwiched between 2 bulk phases Q. 
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However, as already noted in (H [9] [TO ] , the parameters in 
Eq.[4]can be chosen to give a reasonable fit between for 
the relation between grain boundary energy and misori- 
entation for much larger angles than those considered in 
the original problem. 

It is reasonable to assume that the form of Eq. [4] will 
remain valid for binary alloys, with the coefficients E 
and A modified by the presence of segregated solute, as 
well as by the degree of undercooling. This hypothe- 
sis is consistent with phase field and phase field crystal 
simulations of a pure material Ref. (24l[26]|, which used 
Eq.|4]to model grain boundary energy in pure materials 
at different undercoolings by fitting the core energy (i.e. 
parameter A) to temperature. 

With the presence of a grain boundary (or any other 
defect), solute may preferentially segregate to the grain 
boundary. In deriving this modification in the next sec- 
tions, we will be guided by the Gibbs' adsorption the- 
orem, which relates the degree of solute segregation to 
the grain boundary energy and chemical potential of the 
system according to 

where fi x is the chemical potential of species x in the 
binary alloy (x = 1, 2) and T ex = N£ x /A, where N£ x is 
the excess number of atoms of species x in a unit area A 
of the grain boundary (unit length in 2D) Q. It is noted 
that Eq.|5]can be also be written in terms of the chemical 
potential difference, /i^ = fi\ — /i2, if T e x x is replaced by 
TT , the excess particle difference at the grain boundary. 



IV. NUMERICAL RESULTS 

A. PFC Simulations 

We simulate Eqs.[2]and[3]with the Fourier method out- 
lined in Elder and Grant 123) and Mellenthin et al l24l . 
Our computations are performed on a 1024x2048 grid 
with periodic boundary conditions. The grid spacing is 
Ax = 7r/4y/ (1 — 2r]i/;)/(l — Arjip) and the time step is 
At = 1.0. In this work, the average normalized alloy 
concentations studied are ^o=0,-0.05, 0.1, -0.15, and - 
0.2 (-0.2 is only considered for the large angle fits). The 
average normalized density in the PFC model is set to 
no = 0. Bicrystal grain boundaries are formed in al- 
loys whose phase behaviour can be described by a spin- 
odal phase diagram, i.e. w = 0.088, except for some in- 
stances when ipQ = —0.2, for which w = 0.008 | 27 ]. For 
each concentration studied, the parameter Bq is chosen 
such that the undercooling is sufficient for grain bound- 
aries to close (that is, the disjoining pressure does not 
keep the grain boundaries from closing as it could at too 
small undercooling (TIED). The exact value for Bq de- 



pends on the concentration, however, the undercooling 
typically varies from —0.02 to —0.06. The other param- 
eters used in Eqs.[2]and [3] are B% = -1.8, B x = 1, 
t = 0.6, v = 1, u = 4, and K = 4. Note that 77 = 
except where otherwise indicated. 

The basic initial condition for simulating a grain 
boundary begins with two large grains with a small liquid 
gap of roughly 10 grid points in between them. For the 
least deep temperature (Bq ) quenches, at a given average 
concentration ipo, the temperature parameters is dropped 
and a grain boundary forms. Quenches to the lowest val- 
ues of Bq are done in multiple steps, quenching first to 
higher, intermediate temperature, before lowering Bq to 
the final desired value. The crystals are oriented 0/2 and 
— 0/2 from the 0° crystal reference, respectively. Once a 
simulation is started the grains quickly form a bicrystal 
with misorientation, 6. For angles below 4°, the trans- 
lational offset between the grains, which varies for each 
orientation, needs to be chosen carefully to allow for dis- 
locations to form as opposed to the crystals rotating and 
forming one crystal. The grain boundary normal is in 
the x-direction. The angles chosen were: 1.55°, 2.07°, 
2.58°, 3.10°, 4.14°, 5.17°, 6.20°, 8.23°, 10.3°, 12.5°, 
15°, 17.5°, 20°, 22.5°, 25°, 27.5°, 30°, 32.5°, 35°, 40°. 
The small angles are chosen such that an integral num- 
ber of evenly- spaced dislocations fit within the numerical 
domain, so that the results can be compared accurately 
to Eq.|4] Due to periodic boundary conditions, two grain 
boundaries form. The initial spacing between the two 
grain boundaries is 512 grid points, which is chosen such 
that interaction between two grain boundaries is negligi- 
ble (except at possibly 1.55°). Examples of a low angle 
and a high angle grain boundaries are shown in Fig. [T] 

Grain boundary simulations at each Bq were typi- 
cally run for 100000 time steps. Equilibration of a grain 
boundary is determined by the standard deviation of the 
chemical potential, a^ x , and its magnitude relative to the 
mean of the chemical potential, fi x , where x is ip or n. If 
either cr^ c or cr^ n is much greater than 10 -5 or a flc / n c or 
^/x™ / is greater than 10 -2 , simulations are run longer, 
until the above criteria are met. 



B. Calculation of Grain Boundary Energy 

Grain boundary energy is given by the excess grand 
potential per unit area (or length in 2D), 

- (f s - ji^ s - n n n s )) (6) 

where L y is the length of the grain boundary, w g b is the 
width of a region encompassing the grain boundary, f g b 
is the free energy density of the entire domain (includ- 
ing the grain boundary), f s is the free energy density of 
the bulk solid, ifj g b and ip s are the average normalized 
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FIG. 1. Images of concentration (left frames) and density fields 
(right frames) around low angle (bottom frames) and high angle 
(top frames) grain boundaries. The control area enclosing the 
grain boundary lies between the two red lines. The dislocations 
in the low angle grain boundary are shown with yellow squares. 
Note that the x-direction in simulations is the vertical direction 
in the figures. 



concentrations in the grain boundary region and in the 
bulk solid, respectively, while n 9 b and n s are the average 
normalized atomic densities in the grain boundary region 
and in the solid, respectively. i/j g b and ip s and the anal- 
ogous normalized density quantities (obtained by substi- 
tuting i/j with n) are calculated as follows: 



*l>gb 



1 r L y r 

LyW gh Jq J x 

0.5L x ipo 



Xgb+Wgb/2 



ipdxdy 



gb-W gb /2 



0.5L X - w gb 



Given that the ijj s and n s do not differ substantially 
from the average values tpo and no, respectively, f s can 
be written as a Taylor series about no and ^0 as per l24l . 
which gives 

Igb = W gb {fg b ~ ^(l/jgb ~ ^O) 

- fJ>n(n 9 b - n ) - f 8 (n , ^0)} (7) 

where f s (no , ^0) is the f ree energy of a bulk solid at den- 
sity no and concentration ipo, and 
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yWgb Jo J Xgb - Wgb /2 



Xg h +Wg h /2 



fdxdy 



To identify the region encompassing the grain bound- 
ary, the position along the x-axis (transverse to the grain 
boundary) with the maximum average free energy den- 
sity is first found (x g b). Properties evaluated over the 



entire domain, which are needed to determine the extent 
of segregation in the system, are then summed over all 
positions within ± 80 grid points of this reference posi- 
tion (n.b., the effective width of the system is thus taken 
as Wgb = 161Ax, because there are few differences ob- 
served when calculating the degree of segregation with 
Wgb = 141 Ax, 181Ax or 301Ax). 

As noted by 124) . the grain boundary energy is dif- 
ficult to determine because f s is difficult to determine. 
The system parameters were chosen to reduce the effect 
of residual stresses in the bulk on the calculation of the 
solid free energy, although this effect is observed to give 
a relatively uniform contribution to the energy (with the 
exception of when the lattice fits perfectly in the trans- 
verse direction, the bulk solid free energy varies as a ran- 
dom distribution with maximum spread less than 2% of 
its mean). Following (24), we take w g b = L x /2 when 
calculating the grain boundary energy, meaning that Eq 
|7]becomes to first order: 



27^6 = L x {fgb - f s (n ,ipo)} 



(8) 



For each quench temperature Bq and average alloy 
composition ipo, Eq. [4] was substituted for 7^5 in Eq. [8] 
The parameters f s and A were found by fitting the re- 
sulting expression for grain boundary energy versus mis- 
orientation to our computed data. Two types of fitting 
were performed. The first case considered only low 
angles. We assumed a theoretical value for the elas- 
tic modulus Eq analogously to Ref. [24 ] and obtained 
a best fit for f s and A (or alternatively 7*0). The elas- 
tic modulus in this case was taken from the analytic ex- 
pression recently derived in Refs. (191 [23) using a one- 
mode approximation to the total density, i.e. by writing 

n = (j) ^2 j exp (^Gj • x^j , where Gj are the reciprocal 
lattice vectors of the crystal symmetry being considered 
(here 2D HCP). As shown in Ref. 03, *2 for a binary 
alloy is given by 



Y 2 = SB 



Xi2 



(9) 



t+Jt 2 -lbv {AB 



■4BoVWo 2 } 
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where AB Q = Bq — B x and where <p is the equilib- 
rium amplitude of the first modes of the total density n. 
Equation [4] is compared to the computed grain bound- 
ary energies in Fig. [2] only for bi-crystals simulated with 
w = 0.088. The coefficients Eo, A and f s were deter- 
mined as described above. In Fig. [2] the datasets are 
all compressed/stretched by a factor C = ro exp(0.5)/a 
to ensure that all points are fit by the normalized Eq. [4] 
where E = Y 2 b/ (%Tra(j) 2 Wgb) and A = 1.5 - ln(2ir), 
the latter of which is equivalent to ro ~ 4.4. The core 
radii, plotted in Fig. [3] can be seen to be relatively con- 
stant, though possibly decreasing for larger undercool- 
ings Bq — Bq s (Bq s — Bq defines the temperature of 
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FIG. 2. Scaled grain boundary energy, uj g b — 
^gbAxOc/ \wgb4> 2 ), vs. normalized bicrystal crystal mis- 
orientation (in radians for the reference curve), 0/6 c , for 
different concentrations and temperatures for low angles, 
where C = roexp(0.5)/a is a horizontal stretch factor. 
Theoretical reference Read-Shockley curve (solid); = 0, 
Bq = 1.002 (squares); ip = 0, = 0.962 (inverted trian- 
gle); = -0.05, Bq = 0.996 (right triangle); ^ = -0.05, 
Bq = 1.006 (circle); = 0.1, B^ = 1.015 (left triangle); 
^ = 0.1, Bq = 0.995 (star); = -0.15, Bq = 1.035 
(bow tie); Vo = -0.15, Bq = 1.015 (triangle); Vo = -0.15, 
#o = 0-005 (diamond). 



FIG. 4. Scaled grain boundary energy., uj g b = j g bE™(BQ = 
1.002, ipo = 0)/(w gb AxES l (Bfr ,ip )), vs. bicrystal crystal 
misorientation, #/# c (in radians for the reference curve) for dif- 
ferent concentrations and temperatures. Theoretical reference 
Read-Shockley curve (solid). Data for w — 0.088: ipo = 0, 
Bq = 1.002 (squares); ip = 0, Bq = 0.962 (inverted trian- 
gle); Vo = -0.05, Bq = 0.996 (right triangle); ip = -0.05, 
Bq = 1.006 (circle); = 0.1, = 1.015 (left triangle); 
^ = 0.1, Bq = 0.995 (star); Vo = -0.15, = 1.035 
(bow tie); ^ = -0.15, Bq = 1.015 (triangle); ^ = -0.15, 
£^ = 0.995 (diamond),^ = -0.2, Bq = 1.045 (graduated- 
shading circles). Data for w = 0.008: ip = -0.2, Bq = 
1.065 (graduated-shading boxes), = -0.2, Bq = 1.045 
(graduated- shading diamonds),^ = —0.2, Bq = 1.025 (line 
with ties) 
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FIG. 3. Core radius versus undercooling for different concen- 
trations calculated from low angle fit of data in Fig. [2] Square 
(blue) - ipo = 0; Diamond (red) - ipo = —0.05; Upside down 
triangle (yellow) - = 0.1; Triangle (green) - ip = —0.15 



the solidus at a given concentration), consistent with the 
trend found in [ 24 ] for pure materials. 

In the second case, all the raw data of grain boundary 
energy versus misorientation in Fig. [4] are scaled against 
a reference curve given by Eq. [4] with Eq = 1 and 
A = 0.362. To match the simulated grain boundary ener- 
gies to the reference curve, a vertical scaling factor E' = 
E?(B% = 1.002, ^ = 0)(Ax)/E^(B^^j )/w gb is 
applied, where E™ are found empirically by fitting the 
data to Eq.|4] (discussed further below). As in Fig. [2] the 



datapoints are scaled horizontally according to a stretch 
factor, C = exp(0.362 — A). No horizontal shift was 
applied to the data with w = 0.088 for all concentrations 
and temperatures (i.e., A = 0.362). Figure [4] also in- 
cludes data from bi-crystals simulated with w — 0.008, 
which is best fit by a variety of A (0.719, 0.467, 0.411), 
depending on undercooling (Bq = 1.065, 1.045, 1.025, 
respectively). In all cases shown in Fig.|4| 77 = 0. 

The case of w — 0.088 and average concentration 
ipo = —0.15 was also used to analyze the effect the 
degree of mismatch r\ on the thermodynamics of grain 
boundaries. From Fig. [5| we see that the grain bound- 
ary energy is virtually indistinguishable for the non-zero 
r] values studied, except for some misorientations where 
the mismatch strains the structure a fair bit locally (most 
noticeable for B^ = 1.035,7? = 0.05 at 4.13 - 8.28°, 
but also noticeable for 77 = 0.1 at some of the higher 
angles), thereby increasing the free energy of the system. 
For 8.28° it was observed that for larger 77 that the bound- 
ary buckles, thereby relieving some of the stress. The 
7/ 7^ case is observed to give only a small change in 
total interface segregation compared to the 7/ = cases. 
This occurs because r\ does not change the amplitude of 
the density field much, meaning that the Read-Shockley 
pref actor should be approximately the same. This, in 
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FIG. 5. Grain boundary energy at different degrees of mis- 
match. Blue squares, cyan ties - 77 = 0, green triangles, yellow 
triangles - 77 = 0.05, cyan ties, orange shaded balls - 77 = 0.1 



turn, lead to similar grain boundary energies for the two 
cases, as discussed further in Section [TV D| ). 

The data of Figs. [4] and [5] was well-described by the 
Read Shockley relationship using the same A for all 
the undercoolings and concentrations studied, differing 
only between the alloys studied (i.e. w = 0.088 or 
w = 0.008). This implies that for a given alloy system, 
all the curves have a maximum at the same angle. For 
these cases, the low angle behaviour of the data should 
should independently obey Eq.[4]with the same theoreti- 
cal expression of Eq defined through Eq. [9] This implies 
that the effective, empirically-fit, E™ should be related 
to the theoretical one by a factor of around 0.65-0.85, 
which is the approximate ratio of A — In 6 between the 
two fitting schemes for 2° <0 < 12° (empirical method 
typically works because of the scatter in the data at low 
angles). The data of Fig.[6j which show that the theoret- 
ical and empirical values of Eq are related via an affine 
transformation with slope of 0.80 ±0.01 (for the case of 
w = 0.088), is suggestive of this conjecture. The lin- 
ear fit in Fig. [6] does not pass through the origin which 
might suggest that this fit is inappropriate for solids with 
small amplitudes. This consideration is not of concern 
because the relation need only hold for amplitudes larger 
than t/lbv, which is the smallest amplitude a solid can 
have in the PFC model (Eq. 19). 



C. Characterizing Solute Segregation in the PFC Model 

With a semi-analytical characterization for 7^5 in 
hand, we can proceed to derive an expression for the seg- 
regation at grain boundaries as a function of undercool- 
ing and alloy concentration using the Gibbs' Adsoprtion 
Theorem defined by Eq. [5] As mentioned above, seg- 
regation controls several important material properties, 
such as grain boundary melting, second phase formation, 
grain boundary diffusion, creep, etc. 

We begin by noting that in terms of the notation used 



FIG. 6. Effective E™ fit from simulation data vs. low- 
angle theoretical Eq. The line of best fit for the points with 
77 = 0, w = 0.088 has a slope of 0.80 ± 0.01 (squares). The 
points for w = 0.008 are triangles. Their E™ is approxi- 
mately lower by a factors of approximately 1.4, 1.1, 1.05 for 
Bq = 1.065, 1.045, 1.025, which is consistent with their dif- 
ferences in A. The points for 77 / are labeled by the upside- 
down triangles 



in the PFC alloy model, the expression for excess solute 
in Eq.|5]can be recast as 



( d~fgb\ = f<hgb 

= f dlgb dip \ 

dip J TV 



dp \ 



-w gb - 



rex 



(10) 



where it is understood that the formal thermodynamic 
variable T is to be identified by and substituted by the 
reduced temperature variable Bq of the PFC model, and 
Wgb is the width of grain boundary region. Both deriva- 
tives in Eq. flO] can be determined by making use of 
EqHland Eq|9l and by noting that fi^ = df/dip, where 
f (tp , (p(ip) , Bq) is the free energy density derived from 
the PFC model from the single mode approximation of 
the density n. This has been derived in Ref. [19], from 
which it can also be shown that — (w + QB^(p 2 — 
24t]Bq (p 2 )ip + uip 3 . The theoretical quantity /i^ dif- 
fers from the simulated /i^, which also implies that there 
is a systematic uncertainty on any calculated quantities 
which use this value. Using the definition of fi^, we ob- 
tain, 



dtp 



(11) 



(w + (6B%cp - 24rjB x 



■2iP)) + 3uiP 2 +d(P/dip 
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and 



dip d<p dip 
= 2(f) 



(A — 



HO)) 

-4rj 2 )ipO(A 



(12) 



■ HO)) 



^t 2 - 15v(Bfr - B x 
Substituting Eqs.[TT]and [l2]into Eq.fTO] yields 



(B% - 4r] 2 )ip 2 ) 



rex 



1.60 ( 



(w - 



A V BX) 
- Arf)^)6{A 



2ipd(j>/dxp)) 
MO)) 



>/i 2 - 15v(B£ - B x 



4rj 2 )ip 2 ) 



3uip 2 
(13) 



The theoretical expression in Eq. [13] is compared di- 
rectly with numerical simulations of grain boundary seg- 
regation. We follow Cahn's method |2] for determining 
excess solute, using the variable n and ip. This yields an 
expression for the excess concentration in a grain bound- 
ary: 



rex 



^gb 



(l + n) fl 



(14) 



Figure [Tjcompares the segregation predicted by Eq. 13 
to the simulated excess measured via Eq. (14). Eq. ( 13} 
uses the effective Read-Shockley coefficients found in 
the previous section to cover the entire range of sur- 
face energy versus angle. The comparison is very good 
for some system parameters for a broad range of angles 
(e.g., ipo = —0.05, —0.2) and less good, but still approx- 
imately proportional for others (e.g., ipo = —0.15,0.1). 
To estimate the extent discrepancies result from the theo- 
retical analytic one-mode approximation, the relative dif- 
ference between the analytic and simulated fi^ is used to 
determine the relative uncertainties on the analytically- 
determined values. Fig. [5] shows the same compari- 
son, with Eq. 13 re-derived assuming the theoretically 
determined Read-Shockley parameters in the first line 
of Eq. 
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This comparison agrees well with the nu- 
merically determined segregation measurement only for 
small angle data. This is to be expected since the theo- 
retically determined Read-Shockley parameters are only 
rigorously expected to hold at low angles. 



D. Effect of Parameters on Grain Boundary Segregation 

In this section we elucidate the effect of undercool- 
ing, alloy concentration and alloy system on segregation 
to grain boundaries. A typical plot of the effect of vari- 
ous parameters is shown in Fig. [9] Temperature (ABo) 
has a strong effect on the degree of segregation for some 




£ -1.50 -1.00 



r (measured) 



FIG. 7. Prediction of TT from Eq. 



13 



the same quantity 

obtained by direct numerical simulation using Eq. [14] Uses 
empirically determined Read-Shockley parameters. Blue - 
1.015, ipo = -0.15; red - Bq = 1.035, i/j = -0.15; 



-0.15; maroon - Bfi = 1.006, 



B£ 

green - Bq = 0.995, ipo 
^ = -0.05; yellow - Bq = 0.996, ipo = -0.05; cyan - 
Bq = 0.995, ip = 0.1; dark green - Bq = 1.015, ipo = 0.1; 
shaded square - Bq = 1.045, ip = -0.2; shaded ball - 
Bq = 1.015, ip = —0.15,77 = 0.05. A reference line with 
a slope of -1 is included. 




FIG. 8. Theoretical prediction of T^ x from Eq. Ovs. the 
same quantity obtained by direct numerical simulation using 
Eq. [14] Uses theoretical Read-Shockley parameters. Blue - 
Bq = 1.015, ip = -0.15; red - B^ = 1.035, ip = -0.15; 
green - Bq = 0.995, ipo = -0.15; maroon - Bq = 1.006, 
<0 O = -0.05; yellow - Bq = 0.996, ip = -0.05; cyan - 
Bq = 0.995, ipo = 0.1; dark green - Bq = 1.015, <po = 0.1. 
In all cases 77 = and w = 0.088 on this plot. A reference line 
with a slope of -1 is included. 



concentrations, as seen in the difference between the red 
and green curve. This is due to ABo strongly affecting 
the density amplitude, </>, which in turn also affects /i^. 
Furthermore, we see that lattice mismatch has a relatively 
small effect for small \rj\ (< 0.1) as seen in the difference 
between the blue and red curves. However, the effect is 
larger for 77 = 0.2, as seen by comparing the purple and 
red curves. We also see that segregation is much stronger 
for the w = 0.008 alloy system (black curve) than for the 
w = 0.088 alloy system. 
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x 

CD 




FIG. 9. Segregation at misorientation of 22.5° for different 
average concentrations using Eq. 12. Red line- Bq — 0.01; 



Blue circle- rj = 0.1, Bq 



0.2, Bfr = 0.01; Green square- Bq = 



0.01; Purple diamonds- r\ — 
0.01; Black crosses- 



0.01, w = 0.008. 



Physically, we could expect a large size mismatch be- 
tween atoms, r] to alter grain boundary segregation be- 
cause the strains induced in the system by mis-matching 
local lattice constant increases the elastic free energy. To 
reduce free energy, solute atoms typically segregate to 
the grain boundary, which is typically less-ordered than 
the crystal. For smaller mismatch, the effect is not easily 
discernible. The small effect for rj < 0.1 can be ex- 
plained in terms of </> 2 and d\i^ / dip not changing sub- 
stantially with a small change in r\ (these quantities de- 
pend on f] 2 ). This effect can be visualized for low angle 
grain boundaries by considering the Cottrell atmosphere 
around a dislocation core. Assuming linear elasticity is 
valid and that the concentration changes do not strongly 
affect the stresses, the change in concentration far from a 
dislocation decays as l28l : 



AC: 



J xnGbil + v)sin(a) 



1 



(15) 



where r is the distance from a dislocation core, E is the 
Young's Modulus, G is the shear modulus, where a is the 
angle with respect to the slip plane, v is the Poisson ratio 
and X = W> + 1)(1 - il>)/pkT(l + (<91n(a^)/<91n((l + 
ip)/2))), where = exp(— fi^ + /VoVWO is me ac_ 
tivity. In our simulations, 2xt] 2 E is small, so it is ne- 
glected from the expression. The predictions of Eq. 15 
for the solute distribution field were found to approxi- 
mately superimpose on the simulated solute field through 
the core. This simulated solute redistribution can be seen 



FIG. 10. Left: segregation to dislocation core with no mis- 
match, from part of a grain boundary at 4° misorientation for 
^ = -0.05, j] = 0.0, B L = 0.996. Right: Segregation to 
dislocation core and Cottrell Atmosphere around a dislocation 
with mismatch, from part of a grain boundary at 4.14° misori- 
entation for ip = -0.05, r\ = 0.1, B L = 0.996. Blue is small 
ip and yellow is large ip. For reference, there is a bright hollow 
circle to the right of the yellow/red dislocation core in each im- 
age. In the image with misfit, there is a darker blue region in 
the reference circle, indicating solute depletion there. The de- 
pleted solute left of the core was moved around the core, mak- 
ing the yellow/red region bigger than in the simulation without 
mismatch. 



in Fig. 10 for the case r] = 0.1 (right) and r] = (left). 
For larger mismatch, we expect there to be a noticeable 
effect on segregation as shown for r] = 0.2 (purple curve) 
in Fig. [9] This was tested for high angles for ip — —0.05, 
where T^f increased by about 20% with a change of r\ 
from to 0.2, as expected. 

We found that eutectic systems exhibit stronger seg- 
regation than those with a lens-shaped (or double-lens- 
shaped) phase diagram. This is due to the parameter w 
in the binary PFC model, which is related to the differ- 
ence of interspecies bond energy ei2 and self bond ener- 
gies en, 622 ; that is w ~ (2ei2 — en — €22), where 1 
and 2 represent the two kinds of atoms present. As w be- 
comes smaller, the self-attraction becomes larger, mean- 
ing that atoms of type 1 want to be near type 1 atoms and 
type 2 atoms near type 2. This means that a larger de- 
gree of segregation in that this material is more energet- 
ically favourable. This is consistent with what Hondros 
and Seah who related solute enrichment factors (at the 
grain boundary) with solute solubility (which is related 
to 2e i2 - en - e 2 2) 0. 

Another interesting property of the curves in Fig. [9] 
is that a number of the curves have a maximum and a 
minimum degree of segregation (concentration excess). 
We would expect curves in a normal material to have at 
least one extremum because the excess solute is zero in a 
pure material and non-zero value when there is a mixture; 
this does not contradict Hondros and Seah, who only pre- 
dict a monotonic increase in surface excess Fe-P with in- 
crease P concentration because they consider small con- 
centrations (they exclude pure P) so that they can justify 
using Henry's Law to determine the chemical potential 
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l2l . The two extrema (one minimum and one maximum) 
in the PFC model is partially a result of the alloy behav- 
ing mostly like a pure material for ip = 0, so there is 
one extremum between the pseudo-pure (50% composi- 
tion) and each of the pure materials (0% and 100% com- 
position). To predict at which concentration the maxi- 
mum degree of segregation occurs, we could differentiate 
Eq. ( 13 ) with respect to ip. This equation (and a series 
expansion thereof) is rather unmanageable analytically, 
but it is rather easily solved numerically, so this process 
can be adapted to other formalisms. Nevertheless, the 
shift of the peak with temperature observed in Fig. [9] can 
be understood for high angles by considering the grain 
boundary as being an quasi-amorphous state in equilib- 
rium with the solid state. As temperature decreases, the 
amorphous state of the grain boundary becomes more 
metastable relative to the solid, implying that the free en- 
ergy of the grain boundary increases relative to the solid. 
With this relative shift in free energy, a common tangent 
between these two "coexisting" phases steepens, mak- 
ing the intersection with the amorphous phase's "free en- 
ergy" curve (whose minimum is closer to ip = than the 
minimum of the solid) moves closer to zero. 



E. Discussion of Results 

It is instructive to further consider the reasons for the 
weak dependence of grain boundary energy on concen- 
tration in a material where both atomic species are com- 
pletely miscible in each other (i.e. the w = 0.088 
system). As in pure materials, low angle grain bound- 
aries in alloys are expected to have the theoretical Read- 
Shockley-type dependence of grain boundary energy ver- 
sus misorientation. Contributions due to solute segre- 
gation, however, are expected to be minimal; this is 
physically the same as Turnbull's estimation of the in- 
terface energy of a heterogeneous interface, where a 
slight mismatch in lattice parameter leads to little con- 
tribution due to chemical energy (interface between two 
chemically-equivalent phases) |29l . Similarly, we found 
some change in grain boundary energy due to some so- 
lute segregating to the dislocation cores (as indicated by 
Fig. [3]), but because the cores make up only a small frac- 
tion of the volume of the grain boundary, the majority 
of the grain boundary energy would be in the form of 
elastic energy associated with making a line of disloca- 
tions. Since elastic and impurity effect are both mani- 
fested via Eq. [9j other solutal effects on grain boundary 
energy, such as decreasing the core energy, could not be 
separated from the elastic effects accounted for by the 
parameters used in this study. 

On the other hand, the grain boundary energy for 
high angles is expected to be roughly constant Q and 
dominated by the presence of an undercooled metastable 
phase between the crystals [ 30 1. This metastable phase 



is present throughout the entire grain boundary and so 
solute segregation might be expected to show a stronger 
effect for high angle grain boundaries than for low an- 
gle grain boundaries. Fig. [7] confirms that the degree of 
solute adsorption, T^f , increases with angle because the 
larger angles have higher grain boundary energies. How- 
ever, the solute segregation sometimes only decreases 
grain boundary energy slightly. For example, when com- 
paring the various datasets in Fig. |4| we found that for 
the case of w = 0.088, the normalized results for ipo = 
datasets — which are mathematically equivalent to simu- 
lating a pure material — are essentially indistinguishable 
from those of a binary alloy with non-zero ipQ. 

To further demonstrate that the effect of solute segre- 
gation on the internal energy of the w = 0.088 system is 
small, we solved the linearized form of Eq.|3]in a ID ge- 
ometry. The geometry consists of two solid regions with 
an appropriate amplitude for the density oscillations and 
in between them there is a region of uniformly depressed 
amplitude, corresponding to an amorphous phase. The 
widths and amplitudes of the amorphous region are cho- 
sen to give approximately the correct amount of segrega- 
tion. The change in free energy is found to be less than 
10% of the total free energy, which suggests that the ef- 
fect is small. If the degree of solute segregation were 
larger, this approach would have to be applied iteratively 
to both fields before inferring that the effect is small. 

Some experimental and molecular dynamics studies 
have shown that the effect of solute segregation (mod- 
ulated by the material composition) on grain boundary 
energy can be quite large [0 El Q3L contrary to our 
grain boundary energy, which suggests (at least for the 
w = 0.088 alloys) that the change in grain boundary en- 
ergy is mainly due to undercooling (which, of course, is 
a function of ip at constant Bq), rather than solute segre- 
gation directly. Because the systems studied by Kirchner 
and Kieback lfT2l and Seah and Hondros Q have limited 
solubility of solute in the parent phase (i.e. a eutectic al- 
loy), we set the parameter w = 0.008, which simulates 
eutectic phase diagram, to examine if grain boundary en- 
ergy changes significantly between spinodal and eutectic 
alloy systems. Doing so, we observed that the E™ for 
bi-crystals simulated for the w = 0.008 system are lower 
than the those of w = 0.088 system, as shown in Fig. [6] 
This suggests that segregation does, in fact, reduce grain 
boundary energy more significantly in eutectic alloys. 
Repeating calculations with the linearized form of Eq. [3] 
for w = 0.088 further reveals that the degree of seg- 
regation is higher (roughly a factor of 2). Because the 
degree of segregation is much larger for w = 0.008 than 
w = 0.088, solutal effects on the grain boundary energy 
depress 7^5 noticeably — see Fig. [6] 

On the other hand, some materials, where the differ- 
ent atomic species have a strong affinity for each other 
(both species are completely miscible), include alloys 
with lens-shaped phase diagrams (e.g. Fe-Cr at high tern- 
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perature fl3j ) or the double lens-shaped phase diagrams 
(e.g. Cu-Au O). In both these cases, the grain bound- 
ary energy dependence on composition is observed to be 
very weak (if at all present), consistent with the observa- 
tions reported in this work. 

While we only analysed closed grain boundaries in 
the PFC model, Gibbs adsorption theorem, as stated in 
l2l . suggests that it would be applicable even for wet- 
ted grain boundaries. Motivated by the ID PF analyses 
of Mishin et al [6], we checked this behaviour in a ID 
PF model derived from as PFC amplitude model fT9lL 
where the grain boundary is represented as a strained 
region. Using a T, V, ensemble, we can evaluate 
(djgb/dfi^T -I> = - J dxijj - ip s directly by 
varying fi^ (causing p to vary with it), which allowed us 
to verify that a premelted boundary obeys this relation in 
the same way a closed boundary does. This knowledge 
can be a powerful tool for studying premelting in alloys. 

It should also be remarked that while solute segre- 
gation to grain boundaries can have significant effect 
on dendritic growth, where even small changes in sur- 
face energy anisotropy yield measurable changes in mi- 
crostructure morphology. This was recently demon- 
strated by analyzing a PFC model related to the one used 
here in work done by Provatas et. al l22l . Another prop- 
erty that is expected to be more strongly affected by seg- 
regation to grain boundaries and around dislocation cores 
is grain boundary pre-melting behaviour, in which a sys- 
tem might display multiple grain boundary widths at the 
same state variables. These will be investigated using 
PFC simulations in an upcoming work. 



V. CONCLUSION 

The PFC alloy model has been shown in numerous 
works to self-consistently capture the thermodynamics 
and elasto-plasticity inherent in many diffusive phase 
transformations in metals. This work provided another 
test of the robustness of the PFC formalism in predicting 
the important physical property of solute segregation and 
grain boundary energy in binary alloys. 

We used a phase field crystal alloy model to study 
solutal effects on grain boundary properties in spinodal 
and eutectic binary alloys. We derived a semi-empirical 
model of excess solute segregation to the grain bound- 
ary. This model was used to elucidate the role of under- 
cooling, average alloy concentration, lattice misfit on the 
grain boundary energy on spinodal and eutectic alloys. 
We found that for alloys with lens-shaped phase dia- 
grams exhibit a negligible direct contribution to the grain 
boundary energy, both at low and high mis-orientation 
angles. However, undercooling strongly impacts the en- 
ergy in the PFC model through their effect on the elastic 
coefficients. This finding is in agreement with the find- 
ings of other groups, which found that binary materials 



having high solubility of one material in the other show 
little change in grain boundary energy with composition 
changes. On the other hand, we found that for a eutec- 
tic alloy system, solute segregation to the grain boundary 
had a stronger impact on its energy, again consistent with 
other works. 

We also found that (small) lattice mismatches (i.e. Ve- 
gard's law parameter rj) did not strongly affect segrega- 
tion, though we observe that higher degrees of mismatch 
will have a stronger effect on segregation (as predicted by 
13\ and possibly because of that grain boundary energy. 

There are numerous applications of the results and 
methodology found in this article. The phase field crystal 
formalism can be linked to traditional mesoscale phase 
field methods through various coarse graining proce- 
dures [QjO [3lJ [32). As a result, the grain boundary en- 
ergy and segregation results inferred from this work can 
help guide the parameterization of mesoscale continuum 
theories whose forms are often — by necessity — phe- 
nomenological (e.g., (30l|33l). Similar phase field phe- 
nomenologies have been recently used to make predic- 
tions about grain boundary wetting — in particular about 
how the disjoining pressure changes as the system pa- 
rameters change. 

A simpler application of an analysis such as the one 
conducted here might be useful to predict when premelt- 
ing should occur from the grain boundary energies and 
then compare it to the results observed in the binary alloy 
PFC model as was done for a pure material in l24l . Even- 
tually, the PFC formalism should be applicable to more 
alloys with a greater number of constituent elements. We 
hope that this analysis is taken further by systematically 
analyzing the segregation properties and thermodynamic 
effects thereof in many materials; in conjunction with 
theoretical attempts, such as the one in this paper, ex- 
periments need to be done to verify their predictions and 
possibly uncover important, but neglected effects in the 
analyses. 
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